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Abstract 

We study strange dibaryons based on the SU(2)-embedded B = 2 toroidal soliton. Treat- 
ing the excursions of the soliton into strange directions as small rigid oscillations, we obtain 
a good approximation to the bound state approach. We calculate the dibaryon mass for- 
mula to order 1/N and find that the doubly strange I = J = dibaryon is bound by about 
90 MeV. 
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The existence of stable multibaryon states with vanishing hypercharge, although not yet 
experimentally confirmed, remains an exciting possibility. The simplest of such states, the 
doubly strange dibaryon H, was first conjectured to exist in [|TJ. On the basis of an MIT 
bag model calculation, its mass was predicted to be uih = 2150 MeV, well below the AA 
threshold. This means that all strong decays of H are forbidden, and it is expected to have 
a long lifetime typical of weak decays. Experimental detection of doubly strange dibaryons 
is a subtle matter, and there are some remarkable ongoing efforts in this direction H. In the 
meantime it is important to sharpen the theoretical understanding of dibaryons by resorting 
to other available models of low-energy QCD. 

One viable alternative to the quark models is the Skyrme model [f|, |], |5|] where baryons 
are identified with solitons of the non-linear meson action, which is usually taken to be 

r £i i s2 

S = NS WZ + Jd A x ^ti{d^d^U) + ^ti[d^UU\d u UU^f + J -^tiM{U + U^-2) 

(1) 

with U(x,t) e SU(3). N is the number of colors, and the Wess-Zumino term was first 
determined in ||. 

(m\ \ 
M = \ ml 

\ 2m^ -ml J 

is proportional to the quark mass matrix. While m v = 138 MeV is small and is often 
neglected, the effects of = 495 MeV are significant. With the standard fit to the nucleon 
and A masses obtained with this Lagrangian ||, the parameters are assumed to be = 
108 MeV and e = 4.84. 

After a remarkable success in describing the properties of non-strange baryons, based 
on the SU(2) collective coordinate quantization of the Skyrme hedgehog || [7|], there has 
been a number of attempts to model strange dibaryons. The first interesting idea appeared 
in || where an SO (3) imbedded soliton of baryon number B = 2 was constructed. In 
the limit of vanishing m n and mjf, the classical mass of the SO(3) soliton was found to 
be 1658 MeV = 1.92Mh- The classical mass of the Skyrme hedgehog is Mh = 863 MeV. 
Thus, in the chiral limit, the SO(3) soliton is stable against decay into two B = 1 states. 
Its collective coordinate quantization leads to SU(3) multiplets of zero triality and H, the 
SU(3) singlet, is the lightest state ||. An important feature of the SO(3) soliton, however, 
is that it extends significantly into the strange directions of SU(3). As mx is dialed to its 
physical value, which is appreciable, the SO (3) symmetry crucial to the soliton's existence 
is destroyed. A perturbative estimate shows that inclusion of m K pushes M S q(3) well above 
2Mh, destroying its classical stability fll"0"| . Thus, the dibaryon states constructed by the 
collective coordinate quantization of the SO (3) soliton, which were found to be stable in the 
chiral limit, may not survive the breaking of SU(3). This question requires further study. 

In this paper we focus on another B = 2 soliton solution, which is embedded entirely in 
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the light SU(2) subgroup of SU(3), and is based on the cylindrically symmetric ansatz [l I 

/ sin B (r, 9) cos 2c 



iF(r,9)n-T 







n 



sin 6 (r,#) sin 20 j . (2) 
^ cos0(r, 9) 



The classical mass of this soliton has no dependence on mx- By numerically relaxing the 
two unknown functions with the boundary conditions 

F(r = 0) = 7T, F(r = oo) = 0, 6(0 = 0) = 0, 6(0 = tt/2) = tt/2, 



the classical mass was found [llj to be Mb=2 ~ 1660 MeV, which is considerably smaller than 
the mass of the SO(3) soliton evaluated with vtlk = 495 MeV. Constructions of dibaryon 
states based on the soliton (§) may be found in the literature |L2], [T3|, [14]]. In |T^], an SU(3) 



collective coordinate quantization of Ub=2 was carried out, and certain stable dibaryon states 
were predicted. While these predictions are interesting, it should be noted that a similar ap- 
proach to the octet and the decuplet of baryons has not been particularly successful because 
of problems with the breaking of SU(3) [|TIJ. There exists another approach to strangeness 
PI, however, which makes no explicit mention of SU(3) multiplets but nevertheless appears 



to be more successful quantitatively [|I6| , [T7| . In this approach, hyperons are modeled by 
bound states of kaons and SU(2) solitons [[Ljj. 

The bound state approach was first applied to strange dibaryons in ||13|| . To simplify 
calculations, the ansatz (g) was restricted to F(r,9) = F(r), 6(r, 6>) = 9. This restriction 
effectively forces the soliton into a spherically symmetric shape and pushes its classical energy 
above IM^. While an interesting qualitative picture emerged, no definitive assessment of 
dibaryon stability could be made. More recently, an improved study of bound state dibaryons 
was made in ||14|| . In this paper, a better (although not minimum energy) B = 2 soliton 
ansatz was used, and kaon modes were studied in its background. The lightest dibaryon was 
predicted to be bound by about 35 MeV. This calculation, however, did not take full account 
of kaon-kaon interactions, which are difficult to include in the bound state approach, but are 
expected to be particularly important in this system. 

In this letter, we propose a new description of the Skyrme model dibaryons which we call 
the rigid oscillator approach and which is intermediate between the collective coordinate and 
the bound state approaches. As in the collective coordinate quantization, we allow only the 
rigid motions of the soliton. As in the bound state approach, we expand in the deviations of 
the soliton into strange directions, which we denote by D. Since D turns out to be of order 
1 / x/iV, this generates a systematic 1/N expansion. The suppression of the strange deviations 
is related to the fact that, in the quark model language, we are dealing with states consisting 
of 2N — 2 light quarks and only 2 strange quarks. While increasing m# further reduces 
the strange deviations, our methods work for any m^.[] An advantage of the rigid oscillator 

1 This is a Skyrme model manifestation of the general fact that for large N baryons, flavor SU(3) breaking 
is large even for small tuk Pq, E3. 
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approach is its simplicity: it provides analytic formulae for various quantities in terms of 
m K , thus serving as a good physical guide to the bound state calculations. Furthermore, 
it is not hard to include all terms quartic in D, which are necessary for complete order 
1/N calculations. The rigid oscillator approach has been applied to the B = 1 sector in 
fj~8| , p0| , plf , giving a reasonable approximation to the bound state approach. In ||22|| , it 
was further noted that the rigid oscillator approximation gets better with increasing baryon 
density and beyond some critical density becomes identical to the bound state approach. 
We may expect, therefore, that the rigid oscillator approach will work better for the denser 
B = 2 soliton than for the Skyrme hedgehog. 

Let us start by reviewing the rigid oscillator approach to quantization of the B = 1 
Skyrme hedgehog, 

U h (x) - 













While in [|T8| , |20| , pT|j the approximation with = was considered, here we extend it to 
include the pion mass. We consider only the rigid motions of the soliton, 

U(x,t) = A{t)U h {x)A{t) ] . 



To separate the SU(2) rotations from the deviations into strange directions, we write [EC 



A(t) 



A{t) 










S(t) 



where A(t) E SU(2), and 



where 



V 



S(t) = expi d a X a = expiV 

a=4 



V2 






V2D 


V2D ] 






id 5 
id 7 



(3) 



(4) 



(5) 



In calculating to order iV°, we may neglect the dynamics of A(t) and treat the strange 
deviations in the harmonic approximation. This leads to the following effective Lagrangian, 



L = -M h + 4$!^!) + i— (p*D - D j D) - r^ml - ml)D ] D 



(6) 



The quantities $i and Ti, whose integral expressions are given in [2T|, may be evaluated 
numerically, 

$i w 0.00186/MeV , Y x w 0.00398/MeV . 

Canonical quantization of (0) leads to the following Hamiltonian, 



h = M h + -^n f n-i- A 



4$i 



8$i 



N 2 
16$i 



D ] D . 



(7) 
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The order N piece of the Hamiltonian is the classical ground state energy M^. The order 1 
piece includes the terms quadratic in D and IT, and thus may be diagonalized exactly using 
creation and annihilation operators 

where 



m = y/l + lQ{m 2 K - m^r^/N 2 . 

The operators Qr) may be thought of as creation operators for constituent strange quarks 
(anti-quarks). The normal-ordered Hamiltonian to order 1 is given by 

H = M h + u X (Ja + uib ] b ; u x = -^r-(jii - 1) , oil = -^-(^i + 1) . 

Thus, replacing a light quark with a strange quark (anti-quark) costs energy u) x (&i)- Note 
that for rriK = Tn n , u)\ vanishes thereby restoring the original SU(3) symmetry (it costs no 
energy to replace a u or d quark with an s quark) but that ui\ tends to a rather large value, 
N/ (4$i). The Wess-Zumino term, which acts as magnetic field in the D-D^ plane, breaks the 
s «-> s symmetry Using the calculated values of $i and F±, we find that u)\ ~ 200 MeV. This 
is a reasonable estimate of the difference between the strange and the light quark constituent 
masses. It is also a good upper bound on the similar quantity (the lowest mode energy) found 
in the bound state approach to strangeness. We conclude that the rigid oscillator approach 
is a sound approximation and proceed to apply it to the B = 2 soliton. 

Since the B = 2 soliton is less symmetric (its flavor rotations are in general different from 
spatial rotations), we consider rigid rotations both in the flavor space and in the real space, 



U{x,t) = A{t)U B=2 {R{t)x)A{t) 



together with the ansatz of eqs. (§-|5|). To order 1 we obtain an effective Lagrangian very 
similar to that in the B = 1 sector, 

L = -M B=2 + A§ 2 D ] D + iN (d ] D - D ] D) - T 2 (m 2 K - ml)DW . (8) 

The only modifications are the extra factor of 2 in the Wess-Zumino term, and the new 



integral expressions [12 



r_ J 7T 
2 — 

f 2 



§ 2 = Ll ( d 3 x {l - cosF) 
8 J 



} J d 3 x(l - cos F(r, 6)) , 

' f (FF) + (66) sin 2 F + 4 sin 2 F- ^ l " ° 



e 2 f 2 \ r 2 sin 2 9 

where (FF) = (dF/dr) 2 + (dF/d9) 2 jr 2 . Numerical evaluation of the integrals in the B = 2 
soliton background yields 

$ 2 w 0.0038/MeV , T 2 w 0.0079/MeV . 
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Canonical quantization of (H) leads to the Hamiltonian 

h = m b = 2 + ^tfn - iJt- (Dt n - &d) + (r 2 « - ml) + ^ 

Diagonalizing it as before, we find 

N N 

H = M B=2 + u 2 a"a + u) 2 Wb ; ^ 2 = — (/i 2 - 1) , ^2 = tt-(/"2 + 1) : 

4sP 2 4sP 2 

where 



(9) 



fi 2 = ^l+4(m 2 K -ml)r 2 $ 2 /W . 
Substituting the numbers we obtain uj 2 ~ 198 MeV. This is only 2MeV less than U\, the 



corresponding quantity in the B = 1 calculation. In the bound state calculations [13, 14 



the kaon mode energies in the B = 2 and B = 1 backgrounds were also close to each other. 

In order to calculate the 1/N corrections we need to include the terms in the Lagrangian 
which are quartic in D or involve the soliton angular velocities. For B = 1 this was done in 
21 1, and the correction to the Hamiltonian was found to be 

-i- + 2 Cl ?/ • f + c 1 f 2 ) , (10) 
" c=l--^>"l), (11) 



where we omit the subscript 1 throughout the last equation. pf, the isospin relative to the 
body fixed axes, is the momentum conjugate to a which is defined by 

A ] A = -ia ■ t . 
2 

For the non-exotic states containing no 6-quanta, T = ^a^ra. The A-particle has I = P? = 
and T = 1/2, so that the order 1/N correction to its mass is 

Sm A = _£l » 23 MeV , 

where we use fl± = 0.00514/MeV. To order 1/N, our methods give m\ ~ 1086 MeV, which 
is close to its physical value of 1115 MeV. 

In extending the calculation to the B = 2 soliton, in addition to the SU(2) angular velocity 

a we include the spatial angular velocity (3, which is also of order 1/N. The complete order 
1/N correction to the Lagrangian @ is 

1 2 • 1 • • 

SL = -0 2 53(dj + iD ] r,jD - iD^Df + -A(d 3 - 2p 3 + iD ] r 3 D - iD ] t 3 D) 2 

2 3=1 2 

1 ~ 2 2 

+ ~~ x YsiPif ~ NE>] 3-tD- -iN (p ] D - D ] D)D ] D - 2i§2 (D ] a ■ t D - D ] a- f D) 

2 j= i 3 

- (£>+£>) + |$ 2 (Z>+D + D^D) 2 + 2$ 2 (p*D - D^D) 2 + \v 2 {m 2 K - ml)(D^D)' 
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where the additional moments of inertia are numerically found to be 



Vl 2 « 0.0106/MeV , A « 0.0072/MeV , A » 0.016/MeV . 

The isospin and angular momentum relative to the body fixed axes are 

jbj = dL_ jbf = 8L_ 



and we find a constraint 



den 1 Oft 



J? = -2 (l b / + T 3 



The calculation of the Hamiltonian is lengthy, but the result is quite simple, 

SH «» - w 2 + ^ ■ f + + k^' ? m ' • (12) 



where c 2 and c 2 are given by ( O ) with subscripts 2 throughout. Note that, unlike in [|IB|, |TJ 
there is no explicit T 3 2 term in the Hamiltonian. This greatly simplifies the diagonalization 
of the Hamiltonian. 
Using the identities, 

f 2 = Va) 2 + -a t « = -f--0 , 
4 V ; 2 2 V2 J ' 

( JV)2 = j( j + 1 ) ) and (/t/)2 = j(j + i) ) we have 

5M 1/iV = ^- {c 2 K(K + 1) + (1 - c 2 )/(J + 1) + ^p(S 2 - 25)} 

where K = fa + T. The significance of the quantum number K is that K(K+ 1) is equal, up 
to an additive constant, to the quadratic Casimir of the SU(3) representation which emerges 
in the trr = m n limit. In such a limit, c 2 = c 2 = 1 and we recover the mass formula of 
the SU(3) collective coordinate quantization. Furthermore, states with K = merge into 
a (0, N) SU(3) multiplet, states with K = 1 — into (2,N — 1), states with K = 2 — into 
(4, N — 2), etc. 



As discussed in [[L4| , not all possible quantum numbers are allowed. There are certain 
constraints which arise due to special symmetries of the soliton. They work, in effect, to 
insure the correct statistics of the overall wave function. Consider, for instance, dibaryons 
with S = —2 and J = (which immediately implies T = 1 and = 0). Then one 
finds that the allowed (/, K) quantum numbers are (0, 1), (1, 1), (2, 1), (2, 3), etc., while the 
forbidden ones are (1,0), (1,2), (2,2), etc. In the SU(3) limit, the allowed J = multiplets 
are (2,N — 1), (6, iV — 3), etc. Remarkably, these multiplets are also allowed in the quark 
model 0. In general, the states constructed from the B = 2, SU(2) soliton together with the 
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states constructed from the SO (3) soliton appear to cover all the quantum numbers found 
in the quark model. 

Among the S = — 2 dibaryons that we constructed, the lightest one has quantum numbers 
I = J = J3 = and K — 1. Its mass to order \/N is 

M = M B=2 + 2u 2 + ^-^ 2084 MeV . (13) 

The binding energy with respect to AA, calculated entirely within our approach, is (66 + 4 + 
18 = 88) MeV, where we separated the classical, order 1, and order 1/N contributions. The 
lightest S = — 2 dibaryon we found is not the H because it originates from the (2, N — 1) 
SU(3) multiplet, not from the singlet. For arbitrary N it has strangeness S = —2, while 
the H has S = — 2N/3. Thus, for large N the breaking of SU(3) makes the state that we 
found lighter than the H, which appeared in the SO (3) dibaryon quantization. Whether this 
conclusion may be extrapolated to N = 3 is an open question. We believe, however, that the 
quantum numbers predicted by the Skyrme model are in accord with the quark model for 
any N, and that we have presented a good description of strange dibaryons for sufficiently 
large N. 

How well can we trust the binding energy of 88 MeV? Indeed, no model of low energy 
QCD is perfect. In our approach to the dibaryon, quantization of other soliton modes, better 
treatment of the SU(3) breaking, etc., could affect the numbers appreciably. Nevertheless, 
our calculation is attractive for its simplicity, and it gives a good physical picture of both 
B = 1 and B = 2 states. Our results strongly suggest that a tightly bound doubly strange 
dibaryon is theoretically natural. 
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